Role of oxygen-oxygen hopping in the three-band copper-oxide model: quasiparticle 
weight, metal insulator and magnetic phase boundaries, gap values and optical 

conductivity 



O 

(N 



CO 



I 

C 

O 
O 



(N 
> 
O 
&\ 

rn 

ON 

o 
o 



X 



Xin Wang 1 . Luca de' Medici 2 , and Andrew J. Millis 1 
1 Department of Physics, Columbia University, 538 West 120 th Street, New York, New York 10027, USA 
2 Laboratoire de Physique des Solides, Universite Paris-Sud, 
CNRS, UMR 8502, F-91405 Orsay Cedex, France 
(Dated: February 21, 2011) 

We investigate the effect of oxygen-oxygen hopping on the three-band copper-oxide model rel- 
evant to high-T c cuprates, finding that the physics is changed only slightly as the oxygen-oxygen 
hopping is varied. The location of the metal-insulator phase boundary in the plane of interaction 
strength and charge-transfer energy shifts by ~ 0.5 eV or less along the charge-transfer axis, the 
quasiparticle weight has approximately the same magnitude and doping dependence and the qualita- 
tive characteristics of the electron-doped and hole-doped sides of the phase diagram do not change. 
The results confirm the identification of La2CuC>4 as a material with an intermediate correlation 
strength. However, the magnetic phase boundary as well as higher-energy features of the optical 
spectrum are found to depend on the magnitude of the oxygen-oxygen hopping. We compare our 
results to previously published one-band and three-band model calculations. 

PACS numbers: 74.72.-h, 74.25. Gz, 71.10.Fd, 71.27. +a 



I. INTRODUCTION 

The basic structural unit of the copper-oxide high- 
temperature superconductors is the CuC>2 plane and the 
important electronic states are hybridized combinations 
of Cu 3d and O 2p a orbitals. The relevant electronic 
configurations are variations around the Cu 3<i 9 /0 2p 6 
configuration expected in a simple ionic picture. While 
the Cu d 10 configuration is relatively close to the chem- 
ical potential, strong local correlations on the Cu mean 
that the Cu d s state is very far removed in energy, and 
this strongly affects the physics X (Correlations on the 
oxygen have also been argued to be important but in the 
cuprates the probability of having two holes on the same 
oxygen is apparently small enough that these may be 
neglected.) The minimal model encoding this electronic 
structure is the "three-band" modeler— (defined more pre- 
cisely below) which involves the Cu 3d x 2_ y 2 and O 2p a 
orbitals, along with local correlations on the Cu site and 
Cu-0 and 0-0 hybridization. While a qualitative under- 
standing of the model was obtained early on&~— , modern 
theoretical methods, in particular Dynamical Mean Field 
Theory (DMFT)£L11 have added considerably to our un- 
derstanding of the electronic structure of correlated elec- 
tron materials. In particular, single-site DMFT predicts 
that at a carrier concentration of one hole per unit cell 
the three-band model undergoes a paramagnetic metal to 
paramagnetic insulator transition if the charge-transfer 
energy is decreased below a critical value. This critical 
value defines a correlation strength: one may character- 
ize a material as being strongly correlated if the charge- 
transfer energy is such that the paramagnetic insulator 
phase is obtained, and having weak to intermediate corre- 
lation if not. Locating actual materials on this continuum 
of interaction strength is of interest. 

Dynamical mean-field theory has been used by several 



groups to study a simplified version of the three-band 
model, in which the oxygen-oxygen hopping is neglected. 
The pioneering single-site DMFT work of Georges et al*£ 
and Zdlfl et aiJ^ has been followed in more recent years 
by cluster dynamical mean-field studies by Macridin et 
aiji£ and by single-site DMFT studies^ Recent work 
has argued that in modeling specific materials such as 
La2Cu04 it is important to use parameters obtained from 
band theory calculations^— and has stressed, in par- 
ticular, the importance of including oxygen-oxygen hop- 
ping. 

Each of these papers considered one particular set of 
parameter values; however, for a comprehensive under- 
standing of the physics and because it is not clear that 
any one method determines parameters accurately, it is 
important to determine how the physics changes as pa- 
rameters are varied. 

In a previous paper— we presented a comprehensive 
study of a version of the three-band model in which direct 
oxygen-oxygen hopping was neglected. While this study 
provided a qualitatively reasonable account of some as- 
pects of cuprate physics, some aspects of the model were 
found to be in disagreement with experimental data, 
including the particle-hole asymmetry in the magnetic 
phase boundary, the optical absorption strength in the 
region just above the insulating gap, and the dependence 
on do ping of the high energy optical absorption. Refer- 
ences [l7H19l presented results that differed in some re- 
spects from those in our work 2 ^ and argued that the dif- 
ferences arose in part from the use of a more realistic 
band structure, including, in particular, oxygen-oxygen 
hopping. 

In this paper we extend our previous study to consider 
the consequences of oxygen-oxygen hopping. We find 
that inclusion of oxygen-oxygen hopping at values which 
are physically reasonable and consistent with band the- 
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FIG. 1: Sketch of a portion of the CUO2 plane, with the 
coordinate and choice of unit cell in the paramagnetic (solid 
line) and antiferromagnetic (dashed line) cases indicated. 

ory calculations does not appreciably change the metal 
to charge-transfer insulator phase diagram or the princi- 
pal characteristics of the electron- and hole-doped states, 
nor does it resolve the contradiction with data in the 
magnitude of the above-gap conductivity; however, in- 
clusion of oxygen-oxygen hopping does resolve the diffi- 
culties with the magnetic phase boundary and the very 
high energy conductivity. Our results indicate that differ- 
ences between our results and those of Refs. fl7l - [T9l arise 
mainly from differences not in band structure but in cor- 
relation parameters (U and charge-transfer gap) and, in 
some cases, from different results calculated for the same 
parameters. 

The remainder of this paper is organized as follows: 
Sec. |TT] presents the formalism, describing the models we 
studied and methods that we used. In Sec. Mil we present 
numerical results for the phase diagram, staggered mag- 
netization and spectral functions. In Sec. IIVI we present 
the optical conductivities calculated from the three-band 
model and compare them to the conductivities calculated 
for the one-band model. We compare our result to the 
result in Ref. H3 in Sec. [Vj Sec. |VT]is the conclusion. 

II. FORMALISM 

The important structural unit for high-T c supercon- 
ductivity is the Cu02 plane, which in its idealized form 
is a square planar array of Cu02 units. A small portion 
of the plane is sketched in Fig. [TJ We adopt the coordi- 



nate system shown in Fig. Q] and measure lengths in units 
of the nearest neighbour Cu-Cu distance so that the two 
lattice vectors are x and y. Our choice of unit cell in the 
symmetry-unbroken phase is indicated by the Cu02 unit 
enclosed in the solid box. We will also have occasion to 
consider a two sub-lattice antiferromagnetic (AFM) state 
with lattice vectors u = x + y and v = x — y of length 
s/2, for which we choose the unit cell enclosed by the 
dashed box. 

We assume that the electronic physics of the Cu02 
plane is described by the "three-band" model introduced 
in the high-T c context by Emery^ This model involves a 
3d x 2_ y 2 orbital centered on the Cu site and created by 
an operator d\ and two oxygen 2p„ orbitals, p x on the O 
displaced in the x-direction from the Cu and p y on the O 
displaced in the y-direction. The creation operators for 
these states are p xy respectively. In the AFM phase it is 
necessary to distinguish the two sub-lattices. 

The Hamiltonian H = i/band + H mt includes a band 
theoretical part -ffband expressing the possibility of elec- 
tron hopping from d to p orbitals and from p to p orbitals, 
and an interaction term. Following common practice we 
take the interaction part to be a repulsive interaction on 
the Cu d-sites (labeled by j) 

Hint n d-\,jndi,j- (1) 

j 

with U ~ 8— 9eVi 21 ' 22 Because of the low probability that 
any oxygen site has two holes we neglect any interaction 
on the O-sites. 

We write the band theoretic part of the Hamiltonian 
in momentum space. We adopt the Fourier transform 
convention 

= ^E e4fc ' (Hj+p0)c £ ( 2 ) 

k 

with Rj the position of unit cell j, p a the position of 
atom a within the unit cell and N the number of states 
in the crystal. 

For the paramagnetic (PM) phase we order the basis as 
IV') = (dk<T,Px,k<r,Py,kcr), with k a wavevector in the full 
two-dimensional Brillouin zone of the problem and a spin 
indices, and obtain the band theory part of the Hamilto- 
nian as a tight-binding model with a3 x 3 matrix £ i 23 i 24 



H 



PM 



£<1 



2it p d sin ■ 



—2it p d sin -f e p + 2i pp (cos k x — 1) 



-2it 



pd ' 



4t„ 



2it pd sin -f 
4:t pp sin sin 
s p + 2t pp (cosky - 1) 



(3) 



Here t p d is the copper-oxygen hopping, t pp is the oxygen-oxygen hopping. It will be helpful in later dis- 
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cussions to define the p-d level splitting A = e p — Ed, 
which we also vary. 

Different authors have obtained different forms and pa- 
rameters for the downfolded Hamiltonian. For most of 
the paper we use the particular form of oxygen-oxygen 
hopping obtained by O. K. Andersen and collaborators^ 
We take t p d = 1.6 eV and study different values of t pp , 
focusing mostly on t pp = or 1.1 eV— ^ which bracket 
the range of values proposed in the literature! 19 i 25 " — We 
show that most of the low energy physics is not sensitive 
to the precise value of the oxygen-oxygen hopping, while 
the magnetic phase boundary and higher energy optical 
conductivity do depend on this parameter. An alterna- 
tive form for the oxygen-oxygen hopping, with the terms 
proportional to t pp on the diagonal of H being absent, 
has been used by some authors^ ~ 19 ' 25 i 28 ~ — In Sec. IVlwc 
use the Hamiltonian form presented by those authors, 
and the value t p d = 1.41 eV proposed by Refs. Il7l - [l9l 

To treat the Neel AFM phase occurring at and near 
half-filling we divide the lattice into two sub-lattices A 
and B, with the A sub-lattice being that part shown 
within the solid box and the B sub-lattice being the part 
remaining inside the dashed box in Fig. [TJ We Fourier 
transform using Eq. ([5} with appropriate intra-unit cell 
coordinates tied to the magnetic unit cell. In the basis 

W) = (dt,p£ M >P£k*,dL>Px,ka,P*k<r) • thC band the " 

oretic part of the Hamiltonian becomes a 6 x 6 matrix 
that we express in block form as 



jjAFM _ 

The 3x3 matrices are 
Ha = H B 



H A H 

Hjm Hi 



M 



(4) 



turn Monte Carlo impurity solver in its hybridization- 
expansion (CT-HYB) formi 32 i 33 The specifics arc de- 
scribed in Refs. I20II23I . For the phase diagram scan 
(Fig. [2]) and the doping dependence of the quasi- 
particle renormalization factor Z = (1 — dY./dui) -1 
(Fig. IT5j) we used the zero temperature Exact 
Diagonalization(ED) 34 i 35 method. We also cross-checked 
the results of CT-HYB and ED calculations. 

CT-HYB is formulated in imaginary time and to obtain 
real-frequency information we perform analytic continua- 
tion of the imaginary-axis self-energies using the method 
of Ref. H^. From the self-energies we calculate the elec- 
tron Green's function G at frequency z (we choose the 
zero of energy such that the chemical potential fj, = 0) as 



G(z,fc) = [zl-H-£(z,fc)]~ 



(7) 



with X a matrix in which all entries vanish except the 
d-d components. 

The electron spectral function A is 



G(u,k) - Gt(w,fc) 



2i 



(8) 



The optical conductivities are obtained from>2£ 
, n . 2e 2 f°° dLo f d 2 k f{u) - f(u + fi) 

(7(11) ~ 



he a j (2^) 2 n 

x Tr|j(fc)A(w + fi,fc)j(fc)A(w,fe)], (9) 

where c is the c-axis lattice constant, f{oS) is the Fermi 
function, the fe-intcgral is over the magnetic Brillouin 
zone for the AFM case and the full zone in the PM case. 

We will also consider the integrated optical spectral 
weight: 



tpd& 



tr)d.& 2 2t pp COS 



t p de l 2 

£p — ^tpp 
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^2i\jpp cos ' 2 



Ep 2£pp 



^tpp cos k>£ 



-2tr 



-t p de 2 

~ £itipp COS 2" 

COS Ky 



(5) 



(6) 

Later in the paper we will require the current opera- 
tor j(fe) which isa3x3or6x6 matrix in the PM or 
AFM cases respectively. Our Fourier transform conven- 
tion, Eq. (|2j implies that j = 8H/5k x for both PM and 
AFM cases. Tomczak and Bicrmann 3 -^ noted that care 
is required in constructing the current operator: if a dif- 
ferent Fourier transform convention is used, for example, 
omitting the p a terms in Eq. ([2]) , additional terms in the 
current operator must be introduced to recover the cor- 
rect result. These terms are not needed with the choice 
of convention used here. 

We solve the model using the single-site DMF T n i 12 
primarily in conjunction with the continuous-time quan- 



Tic 



a(n')dn' 



(10) 



The fl — > 00 limit, denoted K = K(Q, — > 00), yields the 
familiar /-sum rule)2& 



, 7T 



(2tt) 2 



/(w)Tr 



d 2 H(fc) 

dkl 



(11) 



For comparison we have studied the one-band model 
described (in the AFM phase) by*^ 



ttAFM 

■"■lband 



with 



Ed — 4:t' cos k x cos k y — 2i(cos k x + cos k y ) 
— 2i(cos k x + cos k y ) Ed — 4i' cos k x cos k y 

(12) 



Hint — UcS ^ ndfjUdij- 



(13) 



We take the conventional choice of parameters t = 0.37 
eV, t' = — 0.3i, and choose U c g = 9i to reproduce the 
correlation gap. 
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FIG. 2: Metal-insulator phase diagram in plane of interaction 
U and p-d level splitting A for one hole per Cu02 unit in 
the PM phase at zero temperature, using ED as the impurity 
solver. Solid line (red online): phase boundary calculated 
for t pp = 1.1 eV; dashed line (blue online): phase boundary 
calculated for t„„ = 0.— 






doping 

FIG. 3: Staggered magnetization ((m) = \{n,tt — 'idj.)!) v - s - 
doping in AFM phase. Solid line (red online): t pp = 1.1 eV, 
A = 4.2 eV; dashed line (blue online): t vv = 0, A = 4.5 
eV. Positive doping values indicate electron doping and nega- 
tive doping values indicate hole doping. Common parameters: 
U = 9 eV, t vd = 1.6 eV, and T = 0.1 eV. 



III. RESULTS: PHASE DIAGRAM AND 
ELECTRON SPECTRAL FUNCTION 

A. Paramagnetic phase diagram 

Fig. [5] shows the metal- insulator phase diagram calcu- 
lated at a carrier concentration of one hole per CuC>2 unit 
in the PM phase. There are two important features of 
the phase diagram: A c i , the value of the charge-transfer 
energy at which the insulating solution loses stability and 
A C 2 at which the T = metallic solution loses stability. 
The dashed line (blue online) shows the t pp = case 
studied previously^; the results presented here are gen- 
erally consistent with our previous work but the higher 
accuracy data available to us now have slightly altered 
the phase boundaries. The solid line (red online) shows 
the phase boundary for t pp = 1.1c V, confirming that the 
effect of oxygen-oxygen hopping on the phase boundary 
is small, corresponding to about a 0.5 eV shift in the 
critical A C 2 and a smaller shift in A c i. We note that as 
one turns on AFM order the part of metallic regime that 
immediately follows A c i will have a gap, which we have 
shown in Ref. [2^ to fit the experimentally observed value. 
We focus on this regime in the following discussion. 



B. Magnetization 

Fig. [3] shows the staggered magnetization (m) = 
I (fidrf — ndi) | as a function of doping for the model (solved 
of course allowing for the possibility of antiferromag- 
netism) with t pp = and t pp = 1.1 eV. For each value 
of t pp we chose a A such that the undoped system has 



a gap of 1.8 eV in the AFM phase. The values of A for 
the two tpp are only slightly different: for t pp = 1.1 eV, 
A = 4.2 cV; and for t pp = 0, A = 4.5 eV. This is con- 
sistent with our discussion that t pp shows only a slight 
change on the metal/insulator phase diagram. It is im- 
portant to note that the calculations are performed at a 
relatively high temperature T ps 1200 K; the magnetiza- 
tion is not yet saturated. The dashed line (blue online) is 
the t pp = result: we see a clear particle-hole asymme- 
try, with antifcrromagnetism disappearing at about 0.09 
hole doping and 0.07 electron doping. The solid line (red 
online) shows the t pp — 1.1 cV result: at the electron 
doped side it is almost identical with the t pp = result, 
while at the hole doped side antifcrromagnetism vanishes 
at about 0.06 hole doping. Thus oxygen-oxygen hopping 
has a strong effect on the hole-doped magnetic phase 
boundary. Comparison to experiment is complicated by 
the incommensurate (stripe) magnetism observed on the 
hole-doped side, and its complicated interplay with su- 
perconductivity and the pseudogap. None of these are 
accurately treated by the single-site DMFT approxima- 
tion. 

Our calculations reveal that the magnetic phase 
boundary is affected more strongly by changes in t pp than 
the metal-charge transfer insulator phase boundary. A 
possible reason may be revealed by the observation that 
oxygen-oxygen hopping affects the doping dependence of 
the magnetic phase boundary much more strongly than 
it affects the Neel temperature of the undoped phase. 
As the doping is increased, the calculated Neel tempera- 
ture drops, as does the size of the low-T moment. Thus 
the magnetic transition becomes more weak-coupling-like 
at higher doping, and it is well known that weak cou- 
pling transitions are sensitive to details of the fcrmiology, 
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FIG. 4: Spectral functions for undoped AFM case. A d f and 
Adi shows the majority and minority spin density of states for 
d-orbitals. A p is the total p density of states for a given spin 
(= y4 Px + A Py ) and A to t shows spin averaged total density of 
states [Aot = A p + (A d f + A di )/2]. Upper Panel: t pp = 1.1 
eV, A = 4.2 eV, e d = -8.7 eV, e p = -4.5 eV; Lower panel: 
t pp = 0, A = 4.5 eV, e d = -9.1 eV, e p = -4.6 eV. Common 
parameters: U = 9 eV, t p d = 1.6 eV, T = 0.1 eV. Fermi 
energy is at zero. 

which is in turn sensitive to the oxygen-oxygen hopping, 
whereas the stronger coupling phenomena are governed 
more by local physics, which is less sensitive to details. 



C. Spectral functions 

Fig. U] shows the spectral functions computed for the 
undoped insulating AFM state. We choose parameters 
such that the material is on the metallic side of the metal 
to charge-transfer insulator phase diagram, so that the 
insulating gap arises (in the model) from the presence of 
AFM order. The spectral gap is determined by proximity 
to the phase boundary shown in Fig. [5]Thc small shift in 
the phase boundary with t pp indicated in Fig. [5] would 
lead to a shift in the gap values computed for two differ- 
ent tpp values at the same A and U. Because the value of 
A cannot be determined a priori, we study the effect of 
tpp on the spectral form at fixed physical excitation gap of 
about 1.8 eV, with A slightly adjusted correspondingly. 
With the value of t p d we have chosen here we find that if 
the model has a small enough A to be in the PM insu- 
lating phase then as AFM order is turned on the gap will 
be increased by 0.3 eV or more (the increase is largest at 
tpp = and decreases as t pp increases). 

The bottom panel in Fig. [4] is the t pp = case studied 
in Ref. [2(| We see a gap of about 1.8 eV near the Fermi 
energy. The states just below the gap have an apprecia- 
ble oxygen character and are identified with Zhang-Rice 
singlets while the states above the gap have almost an 
exclusively copper character and are identified with the 
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FIG. 5: Spectral functions for 0.18 hole doped PM case. A d , 
A p and A to t shows respectively spin averaged d, p and total 
density of states. Atot — A v + A d . Upper Panel: t pp — 1.1 
eV, A = 4.2 eV, e d = -7.29 eV, e v = -3.09 eV; Lower panel: 
tpp = 0, A = 4.5 eV, e d = -7.93 eV, e p = -3.43 eV. Common 
parameters: U = 9 eV, t pd = 1.6 eV, T = 0.1 eV. 

upper Hubbard band. Examination of the fe-space struc- 
ture shows that the gap is indirect in all cases, with va- 
lence band maximum at (w/2,iv/2) and conduction band 
maximum at ( 0,7r). T he indirect nature of the gap was 
noted in Refs. ll7H19i The non-bonding oxygen band is 
visible as a delta function centered at w = e p = —4.6 eV. 
The peak at binding energy u> ~ 12 eV represents the 
Cu-d 8 state, which is displaced below the U = 9 eV by 
level repulsion with the oxygen. Turning to the upper 
panel, the t pp = 1.1 eV case, we see that the height of 
the gap-edge peak in the Zhang-Rice region is reduced, 
and the Zhang- Rice absorption is spread out over a wider 
energy range. Also, the structure at higher binding en- 
ergy is altered. The delta function at u> = e p — —4.6 eV 
is replaced by an integrablc singularity and the oxygen 
bands are broadened. The peak representing the Cu-d 8 
state is changed in form, but the t pp evidently has little 
effect on its energy. For smaller U or larger t pp the d 8 
feature may be absorbed into the oxygen bands. 

Calculations of both the hole-doped (shown in Fig. [5]) 
and electron-doped cases (not shown here) show the same 
behaviour that a non-zero t pp expands the non-bonding 
oxygen band to a width of 8t pp , consistent with the above 
discussions. We note that at t pp = 1.1 eV the bands 
change enough to absorb the d 8 peak into the oxygen con- 
tinuum, while if the value oft pp is smaller, e.g. around 0.6 
eV as in Refs. 12514271 . the width of the oxygen band would 
be approximately one half of that of the result of t pp = 1.1 
eV. Fig. [5] shows that the Zhang- Rice singlet is not af- 
fected in the non-zero t pp case, in the sense that the spec- 
tral weights of d-content and p-content remain approxi- 
mately equal. We have also found (not shown, but see the 
discussion of the conductivity below) that the self-energy 
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FIG. 6: Optical conductivity [Panel (a)] and the correspond- 
ing integrated spectral weight [Panel (b)] in the undoped 
AFM case at various t pp values. For t pp = 1.1 eV and 0, 
the parameters are the same as in Fig. [4] For t pp — 0.55 eV, 
the parameters are A = 4.35 eV, Ed = —9.9 eV, e p — —5.55 
eV, U = 9 eV, t pd = 1.6 eV, T = 0.1 eV. 



at low Matsubara frequencies depends on t pp only at the 
5% level, if the distance from the metal-charge transfer 
insulator phase boundary is held constant. Therefore we 
conclude that inclusion of t pp does not have an important 
effect on the characteristics of near-Fcrmi-surface states, 
in particular the Zhang-Rice physics. However it does 
affect the high binding energy features. 



IV. RESULTS: OPTICAL CONDUCTIVITY 

A. Results of the three-band model 

Fig. [5] shows the computed optical conductivity in 
the insulating AFM phase over a wide frequency range, 
along with the integrated spectral weight for three dif- 



ferent values of t pp . (The indirect nature of the gap in 
the single-particle spectrum means that the conductivity 
rises slowly above its onset.) A is adjusted such that the 
gap is around 1.8 eV. In the frequency range between 
1.8 eV and 5 eV the conductivity arises from transitions 
between the upper Hubbard band and the Zhang-Rice 
singlet bands. The detailed line shape depends on t pp , 
in a manner which is roughly consistent with the change 
in spectral function displayed in Fig. 2) The t pp = 
and 0.55 eV cases have a two-peak structure which is 
absent in the t pp = 1.1 eV calculation. The t pp = 0.55 
eV conductivity is sharper than the other two. Some 
of the difference in the traces arise from uncertainties 
in the analytic continuation, but there is a clear trend 
with the integrated spectral weights. As t pp is increased 
the weight in the near gap-edge region decreases, but the 
integrated areas reach similar values by Q, ~ 5 eV [cf. 
Fig-E^b)]. We believe that the robustness of the evalua- 
tion of the gap and similarity of the integrated spectral 
weight means that the main features of our calculations 
are reliable. 

The sharper peaks between 5 eV and 8 eV come from 
the transition between the upper Hubbard band and the 
mainly oxygen band at and below uj = e p . All three 
curves show a two-peak structure in these range. The 
two-peak structure is a reflection of a similar structure 
visible in the upper Hubbard band density of states. This 
structure has been previously discussed in Refs. I36ll40ll4ll 
and the variation in peak height may be traced to the 
changes in the density of occupied states. The integrated 
spectral weights are similar in all curves and we have 
verified that the kinetic energy K(fl) — > oo is consistent 
with the value computed directly from the Matsubara 
axis^ We see that t pp does not have significant effect on 
the conductivity in undoped AFM case. 

Fig. [7] shows an expanded view of the near gap region, 
along with the experimental conductivity^ The magni- 
tude of the conductivity, as noted in Ref. [T7ll20| is still a 
factor of two smaller than experiments^ Ref. Ill states 
that inclusion of apical oxygen bands improves the com- 
parision to experiments; the validity of the Pcicrls phase 
approximation at these energies is also open to question. 
Further investigation of this issue is important. 

Fig. [5] shows the low- frequency part of the optical con- 
ductivity in the undoped AFM case at different A values 
(indicated on the figure), with t pp fixed at 1.1 eV. We see 
that as A is increased (pushing the system further into 
the PM metal region), the gap becomes smaller and the 
magnitude of the conductivity becomes larger; as A goes 
into the PM insulating regime, the gap becomes larger 
and the conductivity magnitude smaller. We see that a 
A between 4.2 eV and 3.7 eV, i.e. slightly smaller than 
A c i but noticeable larger than A C 2 is required to place 
the gap in the experimentally observed range. The cal- 
culations are done at fixed temperature (T = 0.1 eV). 
Because the magnetization is not saturated, the different 
A values imply different staggered magnetizations. In 
general the Neel temperature is lower in the insulating 
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FIG. 7: Low-frequency part of the optical conductivity in the 
undoped AFM case at various t pp values. For t vv = 1.1 eV and 
0, the parameters are the same as in Fig. 2] For t pp = 0.55 eV, 
the parameters are A = 4.35 eV, Ed = —9.9 eV, e p = —5.55 
eV, U = 9 eV, t pd = 1.6 eV, T = 0.1 eV. The experimental 
curve is reproduced from Ref. l42l . 
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FIG. 8: Low-frequency part of the optical conductivity in 
the undoped AFM case at different A values, with t pp fixed 
at l.leV and temperatures indicated. For A = 4.2 eV, the 
parameters are the same as in Fig. [4] For A = 4.7 eV, Ed = 
-8.7 eV, e p = -4.0 eV; For A = 3.7 eV, £ d = -8.7 eV, 
E p — —5.0 eV. Common parameters: U = 9 eV, t pd = 1.6 eV. 
The values of staggered magnetization are 0.59 (A = 4.7 eV, 
T = 0.1 eV), 0.56 (A = 4.2 eV, T = 0.1 eV), 0.47 (A = 3.7 eV, 
T — 0.1 eV) and 0.46 (A = 4.2 eV, T = 1/9 eV) respectively. 



regime than the metallic regime, therefore the magneti- 
zation m for A = 3.7 eV (in = 0.47) is smaller than that 
of A = 4.7 eV (m = 0.59). To make sure that the differ- 
ent magnetization values do not change our conclusion, 
we tune the temperature for A = 4.2 eV calculation such 
that it has similar magnetization as the A = 3.7 eV one 
(shown as dotted line in Fig. |8j). We see that the gap 
essentially does not change, while the magnitude only 



FIG. 9: Comparison of the optical conductivity calculated at 
A = 3.7 eV (between A c i and A c2 ) between AFM and PM 
cases. Parameters: Ed = —8.7 eV, e p = —5.0 eV, (7 = 9 eV, 
t p d = 1.6 cV, tpp = 1.1 eV, T = 0.1 eV. 

changes slightly which may be attributed to the error of 
analytic continuation. 

It is also interesting to consider the effect of antiferro- 
magnetism on the gap magnitude. In studies of a one- 
band model with a bipartite lattice^ a large effect was 
found, with the gap edge shifting up by ~ 40% as the 
PM insulating state was allowed to become AFM. In 
our studies of the three-band model with t pp = 02£ a 
non-negligible, but smaller shift was reported (~ 0.6 eV 
out of 3 eV). Other worbil reported no observable shift. 
Fig. O examines the issue for the case of t pp = 1.1 eV and 
A = 3.7 eV. Antifcrromagnetism does correspond to an 
observable shift of the gap edge to higher frequency, but 
from Fig. [5] one would conclude that the shift is smaller 
(~ 0.2 eV) than in the t pp = case. However, a direct 
quantitative analysis of the gap implied by conductiv- 
ity data such as is shown in Fig. [9] is complicated by 
broadening arising from analytic continuation and high 
temperature. In our previous work>2& we advocated using 
a quasiparticlc analysis, defining the gap edge from zeros 
of dct[cjl — £ — H]. Applying this method to the self- 
energies used in constructing the conductivities shown in 
Fig. yields a gap of 1.5 eV in the PM insulating case 
and 2.0 cV in the AFM case. 

Fig, mil shows the optical conductivities and integrated 
spectral weights for 0.18 hole-doping case. For t pp = 
we had previously found^S that hole-doping (but not 
electron-doping) activated a sharp, large amplitude peak 
at 3.5 eV associated with transitions from the non- 
bonding oxygen band to the Zhang-Rice holes created 
by doping. We see that inclusion of a non-vanishing 
oxygen-oxygen hopping reduces the amplitude of this 
peak (which is not observed experimentally), spread- 
ing the weight over a range of frequencies. Interest- 
ingly we see that the oscillator strength in the "Drudc" 
zero- frequency peak is essentially independent of t pp , 
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FIG. 10: Optical conductivity [Panel (a)] and the correspond- 
ing integrated spectral weight [Panel (b)] in the 0.18 hole 
doped PM case. Parameters for t pp = 1.1 eV: A = 4.2 eV, 
Ed = —7.29 eV, e p = —3.09 eV; Parameters for t pp = 0: 
A = 4.5 eV, Ed = —7.93 eV, e p — —3.43 eV. Common param- 
eters: U = 9 eV, t pd = 1.6 eV, T = 0.1 eV. 



FIG. 11: Comparison of optical conductivities [Panel (a)] and 
integrated spectral weights [Panel (b)] between the three-band 
calculation, one-band calculations and experiments in the un- 
doped case. Solid lines (red online): t pp = 1.1 eV three-band 
calculation, parameters as the upper panel of Fig. [3] Dash- 
dotted lines (magenta online): single-site DMFT calculation 
of the one-band model [Eq. (fl2l ] in the AFM phase. Dashed 
lines (blue online): four-site DC A calculation^ of the one- 
band model [Eq. (|12[l ] in PM phase. Parameters for one-band 
calculations: t = 0.37 eV, t' = — 0.3i, U e s = 9t, Ed = — 4.5t, 
T — O.lt. Dotted lines: experiments of Ref. [421. 



while differences appear at frequencies larger than about 
0.7 eV. This suggests that while a comparison of the 
low-frequency optical conductivity data to model sys- 
tems may be meaningful, the conductivity at higher fre- 
quencies (but still below the charge transfer gap) de- 
pends strongly on details of the model. Note that in 
the single-site approximation employed here the low- 
frequency spectral weight is proportional to the quasipar- 
ticle renormalization Z. The negligible ip P -dependcncc 
seen in the low-frequency part of the conductivity inte- 
gral therefore confirms the negligible t pp dependence of 
Z in this doping range. 



B. Comparison to one- band model 

Fig . ITTI compares the t pp = 1.1 eV conductivity calcula- 
tion of the undoped AFM case to two calculations of the 
conductivity of the one-band Hubbard model [Eq. (|12p]. 
The "1 band" curve (dash-dotted line, magenta online) 
denotes a single-site DMFT calculation of the one-band 
model in AFM phase, with U adjusted to give the same 
gap. One sees that the conductivity at the gap edge 
is much more sharply peaked in the one-band than in 
the three-band model. The "DCA" curve (dashed line, 
blue online) denotes a four-site-cluster DCA calculation 
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of the one-band model in PM phased This calculation 
includes short range AFM correlations, which are seen 
to act further to steepen the line-shape and pile up even 
more weight near the gap edge. Performing a cluster cal- 
culation for the three-band model would be very valuable. 
Despite all these obvious differences, we see that the in- 
tegrated spectral weight of the three cases are similar at 
around 5 eV, above which the oxygen band below e p par- 
ticipates in the optical transition which is not contained 
in any one-band calculations. In other words, an effective 
one-band model describes the degrees of freedom relevant 
to the conductivity up to the scale of ~ 5 eV but even 
in this frequency range the matrix elements and detailed 
structure of the conductivity require information which 
is beyond the scope of the one-band model. 

Fig. [T2] compares the optical conductivity calculated 
from three-band model and one-band model at 0.10 hole 
doping. At this doping value both calculations are in 
PM phase. We see that the one-band model provides a 
reasonable description up to around 3 eV, above which 
transition from the oxygen band to Zhang-Rice band ap- 
pears and the one-band model result deviates from the 
three-band model one. Examination of the insets (which 
are the zoom-in of the low-frequency regime) shows that 
although the detailed line shape of the Drude peaks are 
slightly different, their Drude weights are quite similar. 



V. COMPARISON TO REF. 17 

In this section we turn to a detailed examination of the 
specific parameter values reported in Ref. to describe 
I^-zSr^CuO,}. These parameters place La2Cu04 in the 
strongly correlated, paramagnetic insulator regime of the 
phase diagram. We are interested in the gap value in the 
insulating state, the conductivity in the doped metallic 
state and the doping dependence of the electron quasi- 
particle weight. By comparing these results to experi- 
ment we can assess whether the parameters proposed in 
Ref. [l?] in fact provide a reasonable starting point for a 
description of L^-zSr^CuC^. As a by-product we will 
see that whether one uses the form of the oxygen-oxygen 
Hamiltonian given in our Eq. (|3]) or the form given in 
Refs. IT7MT9IJ25IJ28I - 430I does not affect the physics in any 
significant way. 

We first verify that the model studied here is the same 
as that studied in Ref. [It]. Verification is necessary be- 
cause Ref. [I?] used a tight-binding model which was ob- 
tained from a detailed fit to a band theory calculation and 
differs in two ways from our Eq. ((3]). First, even longer- 
ranged oxygen-oxygen hoppings arc included. However, 
the amplitudes of these long ranged terms are very small 
(less than 5% of the basic t pp ) and (as we shall see be- 
low in Fig. ITB")) have a negligible effect on the oxygen 
density of states. Second, the diagonal elements of the 
hopping do not contain the t pp term so our Eq. ^ must 
be changed so that the diagonal terms become (sd, e p , 
E p ) and the current operator must be correspondingly 
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FIG. 12: Comparison of optical conductivities [Panel (a)] and 
integrated spectral weights [Panel (b)] between the three-band 
calculation and one-band calculation in the 0.10 hole-doped 
case. Solid lines (red online): three-band calculation. Dashed 
lines (blue online): single-site DMFT calculation of the one- 
band model [Eq. (|12[) ]. Both calculations are in the PM phase. 
Insets: zoom-in of the low-frequency part. Parameters for 
three-band calculations: A = 4.2 eV, Ed = —7.64 eV, e„ = 



-3.44 eV, U = 9 eV, t pd = 1.6 eV, t p 



1.1 eV, T = 0.1 



eV. Parameters for one-band calculations: t = 0.37 eV, t' = 
-0.3t, U eB = 9f, e d = -2.42t, T = O.K. 



modified. Ref. [Ttj chose U = 8 eV. The precise parame- 
ter values for t p d and t pp were not given but subsequent 
worbi^ indicates t p d = 1.41 eV and t pp = 0.66 eV. In 
Ref. [13 the values of Ed and e p were obtained by ap- 
plying a double-counting correction to results obtained 
from a band calculation and correspond in our notations 
to A = e p — Ed = 0.34 cV. We have performed calcula- 
tions with the modified Eq. ©, using these parameters, 
at temperatures T = 0.1 eV and 0.05 eV. Comparison 
to Fig. [5] here shows that even for t p d = 1.6, these pa- 
rameters would lie on the insulating side of the phase 
boundary; the smaller t p d would only make the model 
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FIG. 13: Comparison of 0.24 hole doping density of states 
obtained for the model of Ref. [I?] in this work (solid line, 
red on-line) and digitized from Ref. [l?] (dashed line, blue on- 
line), (a) Copper density of states, (b) Oxygen density of 
states. Parameters: t p d = 1.41 eV, t pp = 0.66 eV, U — 8 eV, 
A = 0.34 eV, e d = -5.17 eV, e p = -4.83 eV, T = 0.1 eV 
(this work) and T as 0.008 eV (Ref. [H) . 



more insulating. 

We have performed several tests to verify the consis- 
tency of the models studied here and in Ref. [T3- First 
we have computed the d-level occupancy and total den- 
sity as functions of chemical potential nd(n) and n to t(/^) 
and have verified that they agree with those presented by 
Ref. [13 to within 2%. The projection of our calculated 
many-body density of states onto the oxygen orbitals is 
shown in panel (b) of Fig. [T3l and is seen to correspond 
closely to the results obtained by digitizing the figures 
in Ref. [T3- The Cu-projected densities of states [shown 
in panel (a) of Fig. are qualitatively similar. The 
difference in peak height near lo = is probably a reflec- 
tion of the different temperatures studied. The difference 
in widths of the higher energy features may reflect un- 
avoidable uncertainties in analytic continuation (which 
is insensitive to fine details of high energy features) . We 
note however that the integral from —15 eV to of our 




FIG. 14: Optical conductivity for undoped paramagnetic in- 
sulator calculated for the model of Ref. [fj Parameters: 
Ed = —7.5 eV, e p — —7.16 eV at temperature T = 0.1 eV. 
The arrow indicates the gap value of approximately 2.3 eV 
calculated from the quasiparticle equation of Ref. [361 . 
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FIG. 15: Quasiparticle weight Z = (1 — dT,/dtv)~ 1 \ U j^o as 
a function of doping calculated u sing an ED solver at zero 
temperature, for the model of Ref. fl7l . 



Ad curve reproduces the directly measured d occupancy 
rid ^ 0.53 while the integral of the digitized curve from 
Ref. [13 yields rid ~ 0.58. We therefore suspect that our 
continuation is more accurate. 

We have not been able to stabilize the AFM insulating 
state at the temperatures accessible to us, so we con- 
fine our discussion to PM states. Our calculated conduc- 
tivity for the undoped paramagnetic insulator is shown 
in Fig. [TU We believe that the broadening of the gap 
edge is in part an artifact of analytic continuation, and 
in part a temperature effect. As shown in Ref. |36| the 
gap may be estimated from the difference between the 
highest negative frequency and lowest positive frequency 
at which the quasiparticle equation uj — ReE(w) — = 
is satisfied. This analysis yields a gap of 2.3 eV shown 
as the arrow in Fig. [T3J This gap is larger than the 1.8 
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FIG. 16: Optical conductivities [panel (a)] and integrated 
spectral weight [panel (b)] at 0.24 hole dopin g of two different 
temperatures calculated for the model of Ref. [171. Parameters: 
e d = -5.17 eV at T = 0.1 eV; e d = -5.22 eV at T = 0.05 eV. 
The arrow at the upper right corner of panel (b) indicates the 
kinetic energy K — 0.48 eV independently calculated using 
Eq. (HT|. 



eV found experimentally and it is difficult to see how 
adding antifcrromagnetism to the present calculation will 
do anything but increase the gap. This analysis there- 
fore suggests that the parameters employed in Ref. [l7l 
overestimate the correlation strength in L^CuC^. How- 
ever, it must be emphasized that the extraction of gap 
values from DMFT data in insulating states is compli- 
cated by the intrinsic difficulties of analytic continuation 
(for QMC data) and discretization issues (for ED data) 
and that these difficulties become particularly severe for 
strong correlations (implying large gaps). While we be- 
lieve that the quasiparticle equation method of Ref. [36| 
is the most reliable method for extracting gaps, the issue 
may not yet be settled. 

We next turn to the hole-doped metallic state. Fig. [T5l 
shows the doping dependence of the quasiparticle weight 
Z = (1 - dZ/duj)- 1 computed using the T = ED 
method. The small value and approximately linear dop- 



ing dependence are hallmarks of a strongly correlated 
metal. The temperature dependent ED method (result 
not shown) indicates a weak increase of Z as T is in- 
creased, while analytic continuation of QMC data at 0.24 
hole doping and T = O.leV yields Z ~ 0.16. Fig. 3(b) of 
Ref. [I?] presents results for Z which are inconsistent with 
the results presented here, being larger and less strongly 
doping dependent (Z ~ 0.3 at hole doping x = 0.3, de- 
creasing to Z ~ 0.22 at x = 0.08). The consistency we 
have found, between QMC and ED, leads us to believe 
that the results presented here are correct. 

These results are very similar to those shown at the 
A = 2 eV hole doping curve (black squares) in Fig. 3 of 
Ref. for the model with t pp = 0, indicating that the 
doping dependence is the same (Z ~ x), but the values 
reported for the t pp = case are larger (Z ~ 0.75a:) as 
expected because the t pp = calculation was performed 
for slightly weaker correlations. This indicates that the 
oxygen-oxygen hopping does not have an important effect 
on this aspect of the physics. 

The Z values presented here imply a (zone-diagonal) 
Fermi velocity vp ~ lcV-A at x = 0.24, decreasing ap- 
proximately linearly with doping x. The experimental 
values are around 2eV, independent of x in the range 
0.1 < x < 0.25 and are consistent with values inferred, 
using Fermi liquid theory, from specific heat and quan- 
tum oscillation measurements at dopings x > 0.2^i. 
(Very recent experiments^ find a further 30% decrease 
in vp at low T and uj < 10 mcV; this is still strongly 
inconsistent with our calculations). We believe the dif- 
ference between calculation and data arises because the 
correlations in the actual material are less strong than 
assumed in Ref. UtI 

Panel (a) of Fig. \W\ shows the conductivity calculated 
using analytically continued QMC results at hole doping 
x = 0.24 at two temperatures. Three checks can be made 
to our result. First, we have computed 2dfkr(fi)/7r 
both directly from <r(n) [panel (b) of Fig. [16] and from 
Eq. (TIT]) [finding K = 0.48 eV, shown as an arrow in 
panel (b) of Fig. [T5]. The two results are consistent. 
Second, a straightforward unit conversion indicates that 
our K(u) = 1.8eV) ~ 0.15eV corresponds to N cS ~ 0.28, 
in reasonable agreement with the estimate N e s ~ 0.29 in- 
ferred from Fig. 3(a) of Ref. Il7l Third, the low-frequency 
conductivity is characterized by a "Drude" peak, whose 
area may be obtained (if vertex corrections are negligi- 
ble) from the average of the Fermi velocity over the Fermi 
surface: 

^Drude = J ^\V* F (S)\ (14) 

with s a coordinate along the Fermi surface and v F (s) the 
renormalizcd Fermi velocity (see Appendix for details). 
We evaluated -Korude from the product of the height and 
half-width of the zero frequency peak, obtaining 0.06 eV 
and from Eq. ([P4"|l obtaining 0.07 eV. We have also re- 
peated the conductivity calculation at the lower temper- 
ature T — 0.05 eV; while the form of the conductivity in 
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the mid-IR range somewhat changes, both the area in the 
Drudc peak and the area at f2 < 1 eV remain essentially 
unchanged. 

Although as mentioned above the spectral weight inte- 
grated up to 1.8 eV falls on the curve plotted in Fig. 3(a) 
of Ref. UJY our conductivity differs from that shown in 
Fig. 4 of Ref. [l?], being for example larger by a factor of 
around 2 in the 1-3 eV range. The authors of Ref. [l?] 
inform us 4 ^ that there may be a normalization error in 
the result. 

We now compare the results calculated here to those 
obtained in experiment. From panel (b) of Fig. [TB] we 
see that K(fl = 0.8eV) ~ O.lcV « 1.6A' D rudc- We may 
obtain experimental estimates from Fig. 3 of Ref. |43| and 
from Fig. 4.4 of Ref. H3 (note that -Kband in this reference 
« 0.4 eV). We find K ex P(Cl = 0.8cV) w 0.12±0.02cV and 
^Dmde ~ °- 06 ± O.OleV. We believe that the difference 
between the calculated and measured velocities means 
that the agreement of the conductivity is accidental. 

To conclude this subsection we note that the results 
shown here calculated using the Hamiltonian with a 
slightly different for m 17 i 25 i 28 ~ — than previous sections 
can be compared to the results obtained using the model 
of Ref. d in the rest of this paper. Comparison of Figs. [T3] 
and \TE\ with Figs. [5] and \TU\ reveals that the different 
choice of models does not change the qualitative nature of 
the spectrum and the optical conductivity, and the differ- 
ence can be understood from effectively different values 
of tpp. 



VI. CONCLUSION 

In this paper we have considered the effect of oxygen- 
oxygen hopping t pp on the particle-hole asymmetry and 
optical conductivity of a three-band copper-oxide model 
related to cupratc superconductors. We have found that 
the inclusion of oxygen-oxygen hopping does not change 
the physics in any important way. In particular, in 
the model with large oxygen-oxygen hopping, as in the 
model with no oxygen-oxygen hopping, parameters which 
place the model on the paramagnetic insulator side of 
the metal-insulator phase diagram lead to results for the 
insulating gap and, most importantly, the quasiparticle 
mass, which are in disagreement with experiment, sug- 
gesting that La2Cu04 is not as strongly correlated as may 
have been believed. A model with weaker correlations 
would yield a larger and less strongly doping dependent 
Fermi velocity while vertex corrections (neglected in the 
single-site DMFT but shown via phcnomcnological anal- 
ysis to be important in the cuprate o 37 ' 48 ) could produce 
the correct behaviour of the conductivity. 

Incorporation of oxygen-oxygen hopping does not re- 
solve the most dramatic discrepancy between single-site 
DMFT theory of the three-band model and experiment, 
namely the experimentally observed conductivity near 
the insulating gap edge is much larger than that pre- 
dicted theoretically. Resolving the disagreement prob- 



ably requires the introduction of additional bands, as 
noted in Ref. Il9l . However, inclusion of oxygen-oxygen 
hopping helps resolve a different problem, by broaden- 
ing a peak at ~ 3.5 eV which is theoretically predicted 
to appear upon hole doping. This peak is not observed 
exp er iment ally. 

Inclusion of oxygen-oxygen hopping does have several 
important effects. As previously observed by Ref. [l?], the 
use of a more realistic band structure brings the magnetic 
phase diagram into better agreement with experiment. 
Further, the broadening of the non-bonding oxygen band 
due to oxygen-oxygen hopping may cause the lowest-lying 
d 8 band to be absorbed in the oxygen band rather than 
being visible as a separate feature. 

We have compared our three-band calculations to one- 
band calculations. In the undoped AFM case, the one- 
band model clearly overestimates the conductivity in the 
vicinity of the gap edge; however the optical spectral 
weights in the doped cases at frequencies below ~ 0.8 
eV do agree. This suggests that a reduction of the three- 
band model to a one-band model in AFM case is only 
possible at energies well below the charge transfer gap, 
in agreement with previous remark s 18 ' 49 and in disagree- 
ment with our previous worki 2 ^ In contrast, for PM case 
with doping, the one-band model docs provide reasonable 
description up to around 3 cV. 

We have also compared our results to those presented 
in Ref. [I?]- We find, in disagreement with statements 
in this work, that the main origin of the differences be- 
tween that work and ours is not related to the choice of 
oxygen-oxygen hopping Hamiltonian, but instead relates 
to calculational issues. Our results suggest that the key 
issue in fitting a DMFT calculation to band theory is the 
value of A = e p — e&. The value of A is affected by the 
value chosen for the double counting correction, which is 
not on a firm theoretical basis. 
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Appendix 

In this appendix we present an analysis of the low- 
frequency conductivity, based on the assumptions that 
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the self-energy is momentum independent and the model 
is in a Fermi liquid regime. We shall show that on these 
assumptions the Drude weight is given by an integral of 
the Fermi velocity over the Fermi surface. To do this 
we first note that in the low-frequency limit of the Fermi 
liquid regime the self-energy (in general, a matrix in the 
space of band indices which we denote in bold-faced type) 
can be written 



= — A A 



i-z- 



(A-l) 



Here we have assumed that there is no momentum de- 
pendence, in principle both re-arranges the bands 
and acts as an effective shift of chemical potential; it is 
given by the zero frequency limit of the real part of the 
self-energy. T is the level broadening given by the imag- 
inary part of the self-energy. In Fermi liquid regime T is 
small. Z is a mass renormalization factor which relates to 
the slope of the real part of the self-energy at the Fermi 
energy. 

We may rewrite Eq. ([7]) as 



G _1 (o;, k) = Z- 1 / 2 [u - H*(fc) - »T*] Z" 1 



-1/2 



(A-2) 



with 



H*(fc) = Z 1 / 2 (H(fe) — A M ) Z 1 



/2 



r* = z 1 / 2 rz 1 / 2 

Thus in analogy to Eq. (JSj) we define A* as 

[u - H*(fe) - iT*}' 1 -[uj + H*(fe) + iT*}' 1 



(A-3) 
(A-4) 



A* = 



2i 



= Z- 1 / 2 AZ" 1 / 2 (A-5) 
and the renormalized current operator j* as 

-Z 1 / 2 = - (A-6) 



f(fe) = Z l/2 9H ( fe )^/ 2 - 5H *( fe ) 



dk 



dk 



The the expression of the optical conductivity [Eq. ©] 
can be straightforwardly rewritten as 

, . 2e 2 [°° du f d 2 k f{u) - f(uj + fi) 

a(n) = -ci^-\ W?~ n 

x Tr[j*(fc)A*(w + Q,fc)j*(fc)A*(w,fc)], (A-7) 



Eq. (|A-7[) holds if the self-energy is momentum in- 
dependent and all frequencies are low enough that the 
Fermi liquid approximation is valid. 

Because Z has only positive eigenvalues, the Fermi sur- 
face is the locus of k points k = kp for which H(fe) — A M 
has a zero eigenvalue. Labelling the eigenvalue of H 
which passes through zero at k = kp by we define 
the bare Fermi velocity vp = dsk/dk. To each point on 
the Fermi surface there corresponds a wavefunction ip^ s 
and the Hellman-Feynman theorem implies that 

vp ■ 5k = (Vf |H(fe) - H(k F )\ (A-8) 



with Sk = k — kp. 

Similarly, at the Fermi surface one of the eigenvalues 
E(k) of H*(fe) vanishes and corresponding to this eigen- 
value is an eigenvector 



= z-v^ s /V^i 1 (fe) 



with normalization factor 



z 1 - 1 1 (fc) = (^ s |z- 1 |^ s ) 



(A-9) 



(A-10) 



According to the Hellman-Feynman theorem, if we define 
the renormalized velocity v* F by 



v F -Sk = (^l s \U*(k)-H*(kp)\^l 



FS\ 



then 



v* F (k) 



VF(k) 

Z n \k) 



(A-ll) 



(A-12) 



Note that the renormalization of the velocity with re- 
spect to the band velocity is fc-dependent, even though 
the self-energy is not, essentially because the correlated 
orbital mixes differently with the uncorrelated orbitals 
depending on what Fermi surface point one considers. 

Returning now to Eq. (|A-7I) we see that if T is suffi- 
ciently small then the dominant term in the conductivity 
is obtained by projecting everything onto the Fermi sur- 
face wave function so that (at small 0) 



2e 2 



dcu f d 2 k f(w) - f{w + 0) 



He J_ x 7T J (271)2 n 

x j* s (k)A* FS (io + n, k)ti s (k)A* FS (uj, k), (A-13) 

with j FS ,A FS = (^l s \ j*, A* |*fc S ) the projections onto 
the Fermi surface of the current operator and spectral 
function, "qp" stands for "quasi-particle" . For conve- 
nience we also define T FS = (*£ s | T* 
In the Fermi liquid limit we have 



l*E s > 



j FS (k) = vpjk) 
A FS (uj 7 k) 



FS 



(uj - v* F ■ Sk) 2 + (T* FS ) 2 



(A-14) 
(A-15) 



Note that in Eq. © we assumed that the current is in 
x-dircction thus the x-component of the Fermi velocity 
is taken in Eq. (| A- 14ft . Plugging in Eq. (|A-13j) we have 
(at small 0) 



2T 



FS 



^ + 4(r* s ) 



(A-16) 



Thus just as in the usual case the integral of the low- 
frequency (Drude) part of the conductivity is given by 
the average over the Fermi surface of the renormalized 



14 



Fermi velocity. Writing d 2 k — > kdkdQ we have (at small 
Q) 

, m _2e^ [ d9k F (9)v* F ,M 2 2T* FS 

a ^ l) -hcJ 4^ Km w + 4(r FS f 

= 2e^ r d6k F (6) \v F (9)\ 2T FS 
' he J 4tt 2 2 ^2 + 4 (p* g )2 

The second equality is found by symmetrizing x and y 
directions. Integrating as Eq. (fTO"]) , we find the kinetic 
energy as 

2 f°° fhc\ , , J0 f d6k F {6) . * .... 

Kurude = -J [ ) ^qp(^)^ = J ^2 I ^ ( 61 ) I 

(A-18) 

This can be compared to the Drude weight of the com- 
puted optical conductivity. 

The foregoing is general; in the three-band model of in- 
terest here the matrix Z = diag ((1— while 
A M - diag(-Rc£(w = 0), 0, 0). For the hole doping 0.24 



data shown here, RcS = 3.7 cV. (1 - l^)" 1 = 0.16. 
the Fermi surface implied by this A M has an area corre- 
sponding to a hole doping value of around 0.26. the slight 
deviation from the exact doping value 0.24 is due to the 
temperature/broadening effect. Evaluating Eq. (|A-18() 
yields i^Drude = 0.07 eV. the Drude weight of calculated 
conductivity (shown in Fig. [T6|) is 0.06 eV. The close 
agreement indicates that the low-frequency feature of our 
calculation is reliable. 

We have also done the same calculation for 0.18 hole 
doping, in which case Ed = —5.4 cV, e p = —5.06 eV, 
A M = diag(-3.6eV,0,0), Z = diag(0.11, 1, 1). In this 
case the area enclosed by the Fermi surface indicates a 
hole doping of around 0.06, indicating that at the high 
temperatures that we study, the model is not yet in the 
Fermi liquid regime. We have found that K = 0.04 eV 
from Eq. (|A-18[) , comparing to the Drude weight of cal- 
culated conductivity (not shown) 0.033 cV. Despite the 
deviation from the Luttinger theorem the agreement is 
also good. 
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